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Abstract 

We make use of the recent proof that the critical probability for per- 
colation on random Voronoi tessellations is 1/2 to prove the correspond- 
ing result for random Johnson-Mehl tessellations, as well as for two- 
dimensional slices of higher-dimensional Voronoi tessellations. Surpris- 
ingly, the proof is a little simpler for these more complicated models. 

1 Introduction and results 

The Johnson-Mehl tessellation of R d may be described as follows: particles 
(nucleation centres) arrive at certain times according to a spatial (deterministic 
or random) birth process on M. d . When a particle arrives, it starts to grow a 
'crystal' at a constant rate in all directions. Crystals grow only through 'vacant' 
space not yet occupied by other crystals; they stop growing when they run into 
each other. Also, a new particle that arrives inside an existing crystal never 
forms a crystal at all. This generates a covering of M. d by crystals meeting only 
in their boundaries: every point of M. d belongs to the crystal that first reached it, 
or to the boundaries of two or more such crystals if it is reached simultaneously 
by several crystals. 

These tessellations were introduced by Johnson and Mehl [14] in 1939 as 
spatial models for the growth of crystals in metallic systems. The same growth 
model (but not the resulting tessellation) had been considered earlier by Kol- 
mogorov [15] ; similar models were introduced independently by Avrami [T] , [5] . 
Not surprisingly, these models go under a variety of names (see the references) : 
in mathematics, they tend to be called Johnson-Mehl tessellations, so this is the 
term we shall use here. These models have been used to analyze a great variety 
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of problems from phase transition kinetics to polymers, ecological systems and 
DNA replication (see Evans [8], Fanfoni and Tomellini [9], [10], Ramos, Rikvold 
and Novotny [21], Tomellini, Fanfoni and Volpe [23], [2(5], and Pacchiarotti, 
Fanfoni and Tomellini |21j . to mention only a handful of papers); mathematical 
properties of these tessellations have been studied by Gilbert [T2] , Miles [T7J , 
M0ller [THE], Chiu and Quine [3J, and Penrose [23], among others. 

If all particles arrive at the same time then we get a Voronoi tessellation] 
this was introduced into crystallography by Meijering PjJJ in 1953, although it 
had been studied much earlier by Delesse [5], Dirichlet [7J and Voronoi [27], in 
whose honour it is named. Random Voronoi tessellations have been studied in 
numerous papers: for a list of references, see 0]; here, we shall make heavy 
use of the results in [3j. For a discussion of many aspects of Voronoi and 
related tessellations, including random Voronoi tessellations and Johnson-Mehl 
tessellations, see the book by Okabe, Boots, Sugihara and Chiu [20] , 

In this paper we are mainly interested in random Johnson-Mehl tessellations 
of the plane. As in almost all probabilistic models in the literature, we shall 
assume that the birth process is a time- homogeneous Poisson process of constant 
intensity, say, intensity 1. Thus, our particles arrive randomly on the plane 
at random times t > 0, according to a homogeneous Poisson process V on 
R 2 x [0, oo). 

For each particle arriving at position io£l 2 and time t > we have a point 
z = (w, t) 6 V. The crystal associated to 2= (w, t) reaches a point x G R 2 at 
time d,2(x,w) + t, where d 2 denotes Euclidean distance. (The subscript in the 
notation refers to the power in the norm, not to the dimension: we shall write 
dp for the metric on R d associated to the £ p -norm.) Let || • ||jm denote the norm 
on R 3 defined by 

\\(xi,X 2 ,t)\\j M = ^x\+xl + \t\ = \\{x ll X 2 )\\ 2 + \t\, 

and let d = djM denote the corresponding distance. Then the crystal V z associ- 
ated to z € V may be written as 

V z = { x e R 2 : d((x,0),z) = inf d((x,0),z') }. (1) 

A portion of a Johnson Mchl tessellation is shown in Figure [T] 

In this way we see that the Johnson-Mehl tessellation of R 2 corresponds to 
a two-dimensional slice of the Voronoi tessellation of R 3 defined with respect to 
a slightly unusual metric, c£jm! only the absolute values of the time coordinates 
appear in ([1]), so up to rcscaling (changing the density by a factor of two), it 
makes no difference whether we take the particles to form a Poisson process on 
R 3 or on R 2 x [0,oo). 

In addition to studying Johnson-Mehl tessellations, we shall also study two- 
dimensional slices of the usual random Voronoi tessellation of R 3 : the cells 
are defined exactly as in JT]), but using the usual Euclidean metric d 2 on R 3 . 
In fact, much of what we shall say will apply to more general norms on R 3 ; 
however, the tessellation of R 2 associated to a general norm on R 3 is a rather 
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Figure 1 : Part of a random Johnson-Mehl tessellation of R 2 . The dots are the 
projections onto R 2 of those points z of a Poisson process in R 2 x [0, oo) for 
which the corresponding cell V z is non-empty. 

unnatural object. Indeed, for a general norm, the cells V z C R 2 need not even 
be connected, and the associated graph G-p defined below need not be planar. 
For this reason, we shall focus our attention on tessellations associated to 1 1 • 1 1 jm 
and to the Euclidean norm. Another example we shall consider is the norm l\ 
on R 3 , which may be viewed as l\ ®l\ on R 2 ©R 1 ; the associated tessellation of 
R 2 is a Johnson-Mehl type tessellation in which crystals grow as squares whose 
side-lengths increase at a constant rate. 

As our main focus will be the Johnson-Mehl tessellation, we shall always 
take V to be a Poisson process on R 2 x [0, oo), rather than on R 3 , noting that 
the only effect on the resulting tessellation of R 2 is a rescaling. 

Having defined the cells V z associated to the points z of a Poisson process 
V, there is a natural way to construct an associated graph Gp: the vertex set 
may be taken either to be the set of z S V for which V z ^ 0, or all of V (in 
which case vertices corresponding to empty cells will be isolated) . Two vertices 
are adjacent if the corresponding cells meet, i.e., share one or more boundary 
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points. Ignoring probability zero events, as we may, two vertices are adjacent if 
and only if their cells have a common boundary arc. (Two cells may share more 
than one boundary arc; there is an example in Figure [T]) 

Our aim is to study site percolation on the random graph Gp, or, cquiv- 
alently, 'face percolation' on the tessellation {V z } itself. Let < p < 1 be a 
parameter. We assign a state, open or closed, to each vertex of G-p, so that, 
given V , the states of the vertices are independent, and each is open with prob- 
ability p. We are interested in the question 'for which p does G-p contain an 
infinite connected subgraph all of whose vertices are open?'. Equivalently, we 
may colour the cells V z of the tessellation independently, taking each cell to 
be black with probability p and white otherwise, and we ask for which p there 
is an unbounded black component. We shall switch freely between these two 
viewpoints, writing P p for the (common) associated probability measure. 

Let us say that a point x 6 R 2 is black if it lies in a black cell, and white if it 
lies in a white cell. Note that a point may be both black and white, if it lies in 
the boundary of two cells. Let Co be the set of points of M 2 joined to the origin 
by a black path, i.e., a topological path in R 2 every point of which is black. Let 
z be the a.s. unique point of V in whose cell the origin lies, and let C§ be the 
open cluster of Gp containing z , i.e., the set of all vertices of G-p joined to z 
by a path in the graph G-p in which every vertex is open. Since the cells V z are 
connected (see Section^, the set C C R 2 is precisely the union of the V z for 

z e c G . 

Let 

9(p) = P P (|C G | = oo) = P P (C is unbounded), 

and let 

X (p)=E p (\C°\), 

where E p is the expectation corresponding to P p . Note that the graph G-p 
depends on the metric d as well as on V . Thus 9(p) and x(p) depend on d; most 
of the time, we suppress this dependence. 

We say that our coloured random tessellation percolates if 0(p) > 0. It is 
easy to see (from Kolmogorov's 0-1 law, say) that, in this case, the tessellation 
a.s. contains an unbounded black component, while if 0(p) = 0, then a.s. there 
is none. We write 

Ph = PH(d) = inf{p : 9(p) > 0} 

for the Hammcrslcy critical probability associated to percolation on our random 
tessellation, and 

Pt = Pi(d) = inf{p : \(p) = °°} 

for the corresponding Temperley critical probability. 

Our main aim in this paper is to determine the critical probabilities for 
the Johnson-Mehl tessellation and for a two-dimensional slice of the three- 
dimensional Voronoi tessellation. The corresponding task for the random Voronoi 
tessellation associated to a homogeneous Poisson process on M 2 was accom- 
plished recently in [3j, where it was proved that = Pt = 1/2. 
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Theorem 1. Let d denote either djM or d2, let V be a homogeneous Poisson 
process on M 2 x [0, oo) or on R 3 , and let G-p be the graph associated to the 
tessellation {V z } of R 2 defined by (p}. Then p H (d) = Pr(d) = 1/2. More 
precisely, 9{p) > if and only if p > 1/2 and, for every p < 1/2, there is a 
constant a = a(p) > such that 

P p (size(Co) > n) < exp(— a(p)n) 

for alln > 1, where size(Co) is the area of Co, the diameter of Co, or the number 
\C§\ of cells in C - 

The proof of the corresponding result for Voronoi tessellations in |3J is rather 
lengthy. Much of this proof adapts easily to the Johnson-Mehl setting, includ- 
ing, for example, the analogue of the Russo-Seymour- Welsh Lemma. However, 
the hardest part of the proof, a certain technical lemma, Theorem 6.1 in [3J, 
does not. This result asserts that one can approximate the continuous Poisson 
process V by a suitable discrete process; the proof of this extremely unsurpris- 
ing statement makes up a significant fraction of the length of [5] . The analogue 
of this result for the Johnson-Mehl model is Theorem [5] below; because the 
arguments depend on the details of the geometry, a fresh proof is required here. 
Surprisingly, although the Johnson-Mehl model is more complicated that the 
Voronoi model, the proof turns out to be simpler, though still not short. The key 
difference is that we can use the third dimension of the model to our advantage. 

In the next section we describe basic properties of the Johnson Mchl model. 
In Section [3] we outline the proof of Theorem [TJ assuming Theorem [8l this part 
of the paper consists of a straightforward adaptation of arguments from [3J . The 
heart of the present paper is Section[4l where we prove the technical approxima- 
tion lemma for the Johnson-Mehl model. In the final section we discuss some 
generalizations. 

2 Basic properties 

The probability that some point of the plane is equidistant from four points of 
V is zero. Hence, with probability 1, at most three cells V z of the tessellation 
associated to V meet at any point. We shall always assume that V has this 
property. Similarly, given any measure zero set ./V (for example, the boundary 
of a fixed rectangle), we may assume that no point of N lies in three cells. 
Also, as any ball in R 3 contains only finitely many points of V (a.s. or always, 
depending on the definition of a Poisson process one chooses), we shall assume 
that every disk in R 2 meets finitely many cells V z . 

If we take our metric d to be the Euclidean metric di or the £i-metric d\ (and 
take V to be a Poisson process on R 3 ), then the cells V z are two-dimensional 
sections of (bounded) convex sets (in fact polyhedra) in IR 3 , and hence convex. 
For d = djyi this is not true, but the cell V z associated to a point z = (w, t) G V 
is still a star domain, with centre w: if x £ V z and y is a point on the line 
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Figure 2: A point y on the line segment wx in the plane. As we move towards w 
from x at rate 1, the djM-distance from z decreases at rate 1. The djM-distance 
from any other z' G V decreases at most this fast, so if x lies in V z then so 
does y. 

segment wx, as in Figured then we have 

dju((y,0), z) = \\y — w\\2+t= \\x~w\\ 2 - \\y — x\\2 + t 

= d JM ((%, 0), z) - d JM ((x, 0), (y, 0)) , 

while for any z' G R 3 , 

dju ((y, 0), z') > d 3M ((x, 0), z') - d JM ((x, 0), (y, 0)) 

by the triangle inequality. Since djM((x, 0), z) < gJjm((x, 0), z') for all z' 6 V, 
the same inequality for y follows, i.e., y G V z . Thus V z is a star domain, and in 
particular V z is connected. Of course, the same argument applies to any metric 
rf on R 3 = R 2 © E that is the direct sum of a metric on M 2 and one on R. 

Rather than first constructing a Poisson process V on M 2 x [0, oo), and then 
colouring the points of V black with probability p and white with probability 
1 — p, cquivalently we may start with two independent Poisson processes V + , 
V~ with intensities p and (1 — p), corresponding to the black and white points, 
respectively. This is the viewpoint we shall adopt most of the time. 

In this viewpoint, our state space f2 consists of all pairs (X + , X~) of discrete 
subsets of K 2 x [0, oo). An event E C fl is black-increasing, or simply increasing 
if, whenever Xf ) G E and (X+, A 2 ~) G with A+ C X+ and Xf D X 2 ", 
then (X^~,X^~) G E 1 . In other words, E 1 is increasing if it is preserved by the 
addition of (black) points to X + and the deletion of (white) points from X~ . If 
x G K 2 , then l x is black' is an increasing event, and so is any event of the form 
'there exists a black path Pel 2 with certain properties.' 

It is straightforward to check that Harris's Lemma concerning correlation of 
increasing events extends to the present context; see [3]. 

Lemma 2. Let E\ and E2 be (black-) increasing events, and let < p < 1. 
Then ¥ p (E 1 n E 2 ) > ¥ P {E 1 )¥ P {E 2 ). □ 

Let us note the following simple fact for future reference. 
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Lemma 3. There is an absolute constant A = A(d) with the following property: 
let 5CM 2 be a set with diameter at most s, and let loc(S*) be the event that 
every point of S is within d-distance ^l(logs) 1 / 3 of some point of V . Then 
P(loc(5)) = 1 - o(l) as s -> oo. □ 

This lemma is a simple consequence of the basic properties of Poisson pro- 
cesses (and is also a special case of a very weak form of a result of Penrose [55]); 
we omit the proof. 



3 Reduction to a coupling result 

In this section we present a proof of Theorem [TJ assuming a certain coupling 
result, Theorem [5] below, that will allow us to discretize our Poisson process. In 
a sense, Theorem [5] is a technical lemma, and the arguments in this section are 
the heart of the proof. However, as in [3], the hardest part of the overall proof 
is the proof of Theorem [51 presented in the next section. The arguments in this 
section are, mutatis mutandis, exactly the same as those for random Voronoi 
percolation in [3J, so in places we shall only outline the details. 

Given a rectangle R = [a, b] x [c, d] C R 2 , a < b, c < d, let H(R) = H h (R) 
be the event that there is a piecewise linear path P C R joining the left- and 
right-hand sides of R with every point of P black. When H^R) holds, we say 
that R has a black horizontal crossing. Let V(R) = Vh(R) be the event that 
R has a black vertical crossing, defined similarly. Also, let H W (R) and V W (R) 
denote the events that R has a white horizontal crossing or a white vertical 
crossing, respectively, defined in the obvious way. 

Note that H(R) = 7?b(i?) is a black-increasing event. Also, from the topol- 
ogy of our tessellation, H(R) holds if and only if there is a sequence z\,...,zt of 
black points of V such that the cells V Zl and V Zt meet the left- and right-hand 
sides of R, respectively, and the cells of V Zi and V Zi+1 meet at some point of R 
for each i. 

If no boundary point of R lies in three or more Voronoi cells and no corner 
of R lies in two cells (which we may assume, as this event has probability 1), 
then from the topology of the plane exactly one of the events H^R) and V W (R) 
holds, so F p (H h (R))+¥ p (V w (R)) = 1. 

Note that, from the symmetry of the model with respect to interchanging 
black and white, ¥ p (V vl (R)) = Pi_ p (Vt>(i?)) for any R and any p. Furthermore, 
the metrics d we consider are invariant under rotation (of the plane) through 
tt/2, so ¥ p (H h (S)) = F p (V h (S)) for every square S. It follows that ¥ 1/2 (H(S)) = 
1/2. 

Let f p (p,s) denote the P p -probability of the event H([0,ps] x [0, s]), i.e., 
the probability that a rectangle with aspect ratio p and vertical side length (or 
'scale') s has a black horizontal crossing, and note that 

/ 1/2 (l, S ) = l/2. (2) 

The events H(R) arc defined in terms of the existence of certain black paths 
in a certain black/whitc-colouring of the plane (in which some points are both 
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black and white). This random colouring has the following properties: firstly, 
the event that any point (or given set of points) is black is a black-increasing 
event, so any two such events are positively correlated. Secondly, the distri- 
bution of the random colouring is invariant under the symmetries of Z 2 , i.e., 
under translations (by integer or in fact arbitrary vectors), under reflections in 
the axes, and under rotations through multiples of w/2. Thirdly, well separated 
regions are asymptotically independent: more precisely, let p > and 77 > 
be constants. Given e > 0, if s is large enough, then for Ri and R2 two ps by 
s rectangles separated by a distance of at least rjs, and E\ and E 2 any events 
determined by the colours of the points (of R 2 , not just of V) within Ri and R2 
respectively we have \P(E 1 r\E 2 )-F(Ei)F(E 2 )\ < £• To sec this, note that when 
the event loc(Ri) defined in Lemma [3] holds, the colouring of Ri is determined by 
the positions and colours of the points of V within distance 0((log s) 1 ^ 3 ) = o(s) 
of Ri. 

As noted in [3] , the properties above are all that is needed in the proof of 
Theorem 4.1 of that paper, which thus carries over to the present setting. 

Theorem 4. Let < p < 1 and p > 1 be fixed. 7/liminfs^oo f p (l,s) > 0, then 
limsup s _ f00 / p (p,s) > 0. □ 

Together with ^ , Theorem 2] has the following corollary. 

Corollary 5. Let p > 1 be fixed. There is a constant cq = Co{p) > such that 
for every s$ there is an s > sq with f±/2(p, s) > Co- D 

As in the context of ordinary Voronoi percolation, to prove Theorem [T] it 
suffices to prove the following result, analogous to Theorem 7.1 of [3]. 

Theorem 6. Let p > 1, p > 1/2, c\ < 1 and si be given. There is an s > Si 
such that f p (p, s) > oy. 

Theorcm[T]may be deduced from Thcorcm[6]by using the idea of 1-indepcndcnt 
percolation. The argument is exactly the same as in the Voronoi setting, so we 
shall not give it. 

In the light of the comments above, our task is to deduce Theorem [5] from 
Corollary [5l The basic idea is simple: for s large, we shall show that a small 
increase in p greatly increases f p (p, s) = F p (H(R)), where R is a ps by s rect- 
angle. If H(R) were a symmetric event in a discrete product space, then this 
would be immediate from the sharp-threshold result of Friedgut and Kalai [TT] . 
Unfortunately, H(R) is neither symmetric nor an event in a discrete product 
space, so we have two difficulties to overcome. The first is easily dealt with, by 
working on the torus. 

Let T(s) denote the s by s torus, i.e., the quotient of R 2 by the equivalence 
relation {x, y) ~ (x 1 , y') if x — x' ,y — y' G sZ. Instead of K 2 x [0, 00), we shall 
work in the 'thickened torus' T(s) x [0,t]. Note that we do not wrap around in 
the third direction. It turns out that the precise thickness t is not important in 
the arguments that follow: we could use any thickness t larger than a certain 
constant times (logs) 1 / 3 but bounded by a power of s. For simplicity we shall 
set t = s, working in T(s) x [0, s] throughout. 
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Let us write P ^ for the probability measure associated to a Poisson process 
V on T(s) x [0, s] of intensity 1 in which each point is coloured black with 
probability p and white otherwise, independently of the process and of the 
other points. Alternatively, P^ s ' is the probability measure associated to a pair 
(V + ,V~) of independent Poisson processes on T(s) x [0, s] with intensities p 
and 1 — p, respectively. Our metric d = djM (or d2, or di) induces a metric on 
T(s) x [0, s] C T(s) x R in a natural way. Thus, associated to P^ s ) we have a 
random black/whitc-coloured tessellation of T(s) by the Voronoi cells associated 
to (V, d). 

If we restrict our attention to a region that does not come close to 'wrapping 
around' the torus, then P^ s ) and P p are essentially equivalent. More precisely, 
identifying T(s) x [0, s] with [0, s) 2 x [0, s] C K 3 , we may couple the measures 
Pp^ and P p by realizing our coloured Poisson process on T(s) x [0, s] as a 
subset of that on M 2 x [0, oo). Let e > be fixed, and let R = [es, (1 — e)s] 2 . 
Whenever the event loc(i?) defined in Lemma [5] holds, the colour of every point 
of R is determined by the restriction of the Poisson process to [0, s) 2 x [0, s], so 
the colourings of R associated to the measures P^ s ) and P p coincide. Hence, 
Lemma [3] has the following consequence. 

Lemma 7. Let < a, b < 1 be constant, and let R s be an as by bs rectangle. 
Then for every p we have 

F r M(H(R s )) =F p (H(R s ))+o(l) 

as s — ► co. □ 

Lemma [7] says that when studying crossings of rectangles, we can work on 
the torus instead of in the plane. On the torus, there is a natural way to convert 
H(R) into a symmetric event; we shall return to this shortly. 

As in [5] , we wish to apply a Friedgut-Kalai sharp-threshold result from [TT] . 
A key step is to approximate our Poisson process V on T(s) x [0, s] by a discrete 
process. Given 5 = 8(s) > with s/6 an integer, partition T(s) x [0, s] into (s/5) 3 
cubes Qi of side-length 5 in the natural way. (We may ignore the boundaries 
of the cubes, since the probability that V contains a point in any of these 
boundaries is 0.) As in [3], the crude state of a cube Qi is bad if Qi contains one 
or more points of V~ , neutral if Qi contains no points of V~ U V + = V, and 
good if Qi contains one or more points of V + but no points of V~ . Let 6 = S(s) 
be a function of s that tends to as s — > oo (later, S(s) will be a small negative 
power of s); all asymptotic notation refers to the s — > oo limit. Writing 7 = S 3 , 
since S(s) — > 0, each Qi is bad, neutral or good with respective probabilities 

Pbad = 1 - CXp (-7(1 -pj) - 7(1 -p), 
Pncutral = 63tp(— 7), (3) 

Pgood = exp (-7(1 - p)) (1 - exp(-jp)) ~ 7p. 
Also, the crude states of the Qi are independent. 
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Writing N = (s/8) 3 for the number of cubes, and representing bad, neutral 
and good states by —1, and 1 respectively, the measure P^ s ) induces a product 
measure on the set il n = {— 1)0, 1}^ of crude states. 

To prove results about the continuous process, we shall pass to the discrete 
setting and then back; starting from a realization (P^,Pi) of our Poisson 
process, first we generate the corresponding crude states, and then we return to 
a possibly different realization (P21P2) consistent with the same crude states. 
An event such as H(R) need not survive these transitions: a point x or path P 
may be black with respect to {P\ ,Pi) but not with respect to (P^P^)- To 
deal with this problem, we consider a 'robust' version of the event that a point 
or path is black. 

Given r\ > 0, let us say that a point x £ T(s) is rj-robustly black with respect 
to (V + , P~) if the closest point of P + to x is at least a distance 77 closer than the 
closest point of P~ , where all distances are measured in the metric d. [Note that 
whenever x is rj- robustly black, the entire (r7/2)-neighbourhood of x is black. 
There is no reverse implication: for any 77 > and any r, it is possible for the 
r-neighbourhood of x to be black without x being ry-robustly black.] A path P 
is rj-robustly black if every point of P is 77-robustly black. Set 

Cd = swp{d(x, y) : x,y £ [0, l] 3 } < 00. 

If (P^,P^), i = 1,2, are realizations of our Poisson process on T(s) x [0, s] 
consistent with the same crude state, and a point x £ T(s) is (2Cd5)-robustly 
black with respect to (P^ , P^ ) , then it is easy to check that x is black with 
respect to (P^ , ) . 

Our starting point for the proof of Theorem [5] is Corollary O which gives us 
(with reasonable probability) a certain black path. Fortunately, we can 'bump 
up' a black path to a robustly black path at the cost of increasing p slightly, 
using the following analogue of Theorem 6.1 of [3]. Here, and in what follows, 
we say that an event holds with high probability, or whp, if it has probability 
1 — o(l) as s — > 00 with any other parameters fixed. 

Theorem 8. Let d denote either djM or c?2, and let < p\ < P2 < 1 and 
e > be given. Let 5 = S(s) be any function with < <5(s) < s~ e . We may 
construct in the same probability space Poisson processes , "P-f , P^ and P^ 
onT(s) x [0, s] of intensities pi , 1— pi, pi and I—P2, respectively, so ihatPf 
and P^ are independent for i — 1,2, and the following global event E g \ holds 
whp as s — > 00: for every piecewise-linear path P\ C T(s) which is black with 
respect to (P^,Pi) there is a piecewise-linear path P2 C T(s) which is (2CdS)- 
robustly black with respect to (p2ip2~)> suc ^ that every point of P2 is within 
distance log s of some point of Pi and vice versa. 

The proof of this result is a little involved, and will be given in the next 
section. This is the only part of the present paper that is essentially different 
from the arguments for usual Voronoi tessellations given in [3] . 

Before turning to the proof of Theorem let us outline how Theorem [6] fol- 
lows. The proof is exactly the same as that in Section 7 of [3], mutatis mutandis: 
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we use Corollary [5] and Theorem [5] in place of their analogues Corollary 4.2 and 
Theorem 6.1 of [3]; we write (2CdC>)-robustly black in place of 4<5-robustly black; 
7 = 5 3 , the volume of each small cube Qi, replaces 7 — 5 2 , the volume of a 
small square Si in [3J; finally, A" = (s/5) 3 , the number of cubes Qi, replaces 
N = {s/S) 2 , the number of squares S%. 

Very roughly, the strategy of the proof is as follows (for details see [3J). Fix 
p > 1/2, let e > be chosen below (depending on p), and let s be 'sufficiently 
large'. From Corollary[5l after increasing s, if necessary, the crossing probability 
A/2(10, s/13) is at least some positive absolute constant. By Lemma[71 it follows 
that in the torus, i.e., in the measure pT^i ^ nG probability that a given 10s/13 
by s/13 rectangle has a black horizontal crossing is also at least a positive 
constant. Set p' = (p + l/2)/2, say, so that 1/2 < p' < p, and set 5 = s~ s , 
decreasing 6 slightly if necessary so that s/5 is an integer. Using Theorem [51 
we can convert a black path in P-jyg ^° a ' ncar by' robustly black path in F^,^: 

it follows that the P^, -probability that a given 3s/4 by s/12 rectangle has a 
(2C<i#)-robustly black horizontal crossing is not too small, i.e., is at least some 
constant c > 0. 

The sharp-threshold result we shall need is Theorem 2.2 of [3], a simple 
modification of a result of Friedgut and Kalai, Theorem 3.2 of [IT] . Consider 
the state space f2jy = { — 1, 0, 1}^ with a product measure, in which the coordi- 
nates are independent and identically distributed. An event E in this space is 
increasing if u> = (u>i)f =1 £ E and Ui < lo[ for every i imply lu' £ E. Also, E is 
symmetric if there is a group acting transitively on the coordinates 1, 2, . . . , N 
whose induced action on f2jv preserves E. Roughly speaking, Theorem 2.2 of [3] 
says that if E is a symmetric increasing event in f2jv , and we consider product 
measures on 0,^ in which the probability that a given coordinate is non-zero is 
'small', say bounded by p ma , x , then increasing the probability that each coordi- 
nate is 1 by at least A and decreasing the probability that it is —1 by at least 
A is enough to increase the probability of E from 77 to 1 — 77, where 

A = Clog(l/r/)p max log(l/p max )/logA/' 

and C is constant. (For details, see (5J.) 

To apply the result above, we need a symmetric event in a discrete product 
space. To achieve symmetry, following the notation in Section 7 of [3] , we simply 
consider the event E 3 that some 3s/4 by s/12 rectangle in T(s) has a robustly 
black horizontal crossing. To convert to a discrete product space, we divide 
T(s) x [0, s] into A^ = (s/5) 3 cubes of volume 7 = <5 3 , and consider the crude 
state of each cube as defined above. Let E" udc be the event that the crude 
states of the cubes are consistent with £3, which may be naturally identified 
with an event in Q N . Note that F^ s) (E^ udc ) > F^ s) (E 3 ) > c> 0. 

The Friedgut-Kalai result implies that a small increase in the probability 
of black points increases the probability of £ , g rudc dramatically (details below) . 
It follows that Py (s) (£f udo ) is very close to 1 if s is large. Hence, there is 
a very high Pp'^-probability that some 3s/4 by s/12 rectangle has a black 
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crossing. (Not necessarily a robustly black crossing: in passing to the discrete 
approximation and back again, points of our Poisson process may move slightly. 
However, as noted above, any crossing that was robustly black remains black.) 
By a simple application of the square-root trick, one can deduce that the P^W- 
probability that a fixed s/2 by s/6 rectangle in T(s) has a black horizontal 
crossing is also very close to 1 . Finally, using Lemma [7] again it follows that 
/ p (3, s/6) can be made arbitrarily close to 1, and Theorem [5] follows. 

Turning to the quantitative application of the sharp-threshold result, we may 
take 77 to be a (very small) absolute constant. Each cube Qi is very small, and 
the probability that a cube Qi is either good or bad is at most 7 = <5 3 , so we may 
take p m ax = 7- Also, passing from Pp, (s) to Pp (s) increases the probability that 
a given cube is good, and decreases the probability that it is bad, by roughly 
(jp-p')T, see ©. Hence, to deduce Theorem [5] from the Friedgut-Kalai result, 
we need 

(p-p')7>C" 7 log(l/7)/logiV, 

for some constant C . With p and p' fixed, this reduces to C" log(l/7)/ log N < 
1. Since N = s ^ 1 ' and 7 = s~ e ( £ ', this condition can be met by choosing e 
sufficiently small. Note that it is irrelevant whether 7 = 6 3 and TV = (s/5) 3 , as 
here, or 7 = S 2 and N = (s/5) 2 , as in [3]. Indeed, this part of the argument 
works unchanged in any dimension: the key point is that we can afford only 
to discretize to a scale S given by an arbitrarily small negative power of s. 
Fortunately, Theorem [8] applies for such a S. 

4 Replacing black paths by robustly black paths 

It remains only to prove Theorem [H] Roughly speaking, this states that a small 
increase in the probability p that each point is black allows any black path to 
be replaced by a nearby robustly black path. The proof, to which this section 
is devoted, turns out to be the hardest part of the paper. 

We shall use the following fact about random Voronoi tessellations in three 
dimensions; here O denotes the origin. 

Theorem 9. Let V be a homogeneous Poisson process on K 3 of intensity 1, and 
let d denote di or c?jm- For A > let Ek = Ek.A be the event that V contains k 
points P±, . . . , Pfc with the following property: there are points Qi, . . . , Qk £ K 3 
and real numbers < ri,. < Ak 1 / 3 such that d(0,Qi) = d(Pi,Qi) = r,- L 

for every i, and d(Pj, Qi) > r, for all i and j. If A > and C > are constant, 
then 

¥(E k ) = o(e- Ck ) 

as k — > 00. 

Theorem [5] is essentially equivalent to the following statement: if Vb is the 
cell of the origin in the Voronoi tessellation of M 3 defined using the point set 
V U {0} and metric d, then the probability that Vq has at least k faces is 
o(e~ ck ), for any constant C. In fact, Theorem [5] easily implies this: suppose 



12 



that Vq has fc faces, let Pi, ... , Pk be the centres of the adjacent Voronoi cells 
Vi,...,Vk, and let Qi be a common point of Vo and Vi for each i. Then O 
and Pi are the two closest points of V U {0} to Qi, so cither holds, or wc 
have d(0, Qi) = d(Pi, Qi) > Ak 1 / 3 for some i. But in the latter case there is a 
ball B (defined with respect to d) of radius r = Ak 1 / 3 meeting the origin and 
containing no points of V. Placing O(l) balls Bi of radius r/2 so that any such 
ball B contains one of the Bi, the probability that some Bi contains no points of 
V is 0(cxp(— vol(-B,))J , which is much smaller than e~ ck if we choose A large 
enough. 

Let us remark that the statement above may well be known, at least for 
d = d-j. In two dimensions, very precise results are known; Hilhorst [13] has 
shown that the probability pk that the cell of the origin has exactly k faces 
satisfies 

as k — > oo, where A is a certain constant given as an infinite product. It is very 
likely that similar (c/fc) fe asymptotics hold for Johnson-Mehl tessellations; to 
keep the proof simple, we prove only the much weaker o(e~ ck ) bound. 

Proof of Theorem^ The strategy of the proof is as follows: we shall show that 
if Ef. holds, and certain 'bad' events B\, . . . , B 4 of probability o(e~ ck ) do not 
hold, then an impossible situation arises. In the proof, by 'distance' we mean 
the usual Euclidean distance. Constants in O(-) notation may depend on A and 
C unless explicitly stated otherwise. 

For a > constant, let N a be the number of points z of V with d(0, z) < 
ak 1 / 3 . Then E(N a ) = Q(a 3 k), so, choosing a small enough, E(N a ) < k/8. 
Fixing such an a from now on, let B\ be the 'bad' event 

Bi = {N a > fc/4}. 

The probability that a Poisson random variable with mean fc/8 exceeds is 
fc -e(fc) = o( e - ck ), so V(Bi) = o(e- Ck ). 

Let S be a (Euclidean) ball of radius 9(A: 1 / 3 ) containing {z : d(0,z) < 
2Ak 1 / 3 }, noting that if Ek holds, then all Pi lie in S. Let N be the number of 
points of V in S, so has a Poisson distribution with mean vol(5) = 0(k). 
Setting d = 2vol(5)/fc, let B 2 be the event 

B 2 = {N > Cik}. 

Then V(B 2 ) = fc- e(fe) = o(e- Ck ). 

Given N = \V n S\, we may generate P PI S as a sequence Zi, . . . , zjy, where 
the Zi are independent and each Zi is chosen uniformly from S. Let c be a small 
constant to be chosen below, and let B% be the event that there are at least fc/4 
points Zi G V l~l S with d 2 (zi, zj) < c for some j < i: 

B 3 = { \{ Zl :3j<i with d 2 (zi, Zj) < c}\ > fc/4 }. 
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Given N and Z\, . . . , Zi—i, the probability that d 2 (zi, Zj) < c for some j < i is 
at most pi = (i — l)47rc 3 /(3vol(S')). Conditional on N, the number of points Z{ 
within distance c of an earlier zj is thus dominated by a binomial distribution 
with parameters N and pn, whose mean is at most A r2 47rc 3 /(3vol(S')) = 
<d(N 2 c 3 /k). Whenever B 2 does not hold, we have N < C\k, and so ^jy = 
0{c 3 k). Hence, P(i?2 H P3) is at most the probability that a certain Binomial 
distribution with mean 0(c 3 k) exceeds fc/4. Choosing c small enough, this 
probability is o(e~ ck ), so P(P 3 ) < P(P 2 ) + P(B£ n P 3 ) = o( e - CTc ). 

Let ci be a small constant to be chosen below, and let N' be the number of 
points P 6 V n S for which the angle between OP and the plane is within c\ 
of 7r/4. Then AT' has a Poisson distribution with mean 0(civol(5)) = 0(cifc). 
Choosing ci small enough, E(iV') < fc/8. Let B4 be the event 

Bi = {N' > fc/4}, 

noting that P(B 4 ) = o(e _cfc ). 

Suppose that Ef. holds and that none of the events Bj, 1 < i < 4, holds. 
To complete the proof of Theorem [£l it suffices to deduce a contradiction. As 
Bi, B3 and B4 do not hold, there is a subset X of {Pi,...,Pfe} of size at 
least fc — 3fc/4 = fc/4 with the following properties: for every Pi € X we have 
d(0,Pi) > ak 1 ' 3 , the angle Qi between OPj and the plane is at most 7r/4 — C\ 
or at least n/4 + ci, and no two Pj, Pj £ X are within distance c. 

Projecting the 8(fc) points Pi G X from the origin onto the unit sphere 
centred at the origin, we find two points Pj, Pj 6 A" whose projections are 
within distance 0(fc -1 / 2 ). Without loss of generality we may assume that i = 1, 
j = 2 and that P2 is at least as far from O as Pi is. Let P[ be the point of the 
line segment OP 2 at Euclidean distance d,2(0,Pi) from O; see Figure [31 Then 
d 2 (Pi,P{) = OCfc-Vada^Pi)) = Oik^k 1 / 3 ) = o(l). As d 2 (P,P,) > c, 
for fc sufficiently large it follows that d,2(P{,P2) > c/2, i.e., that d2(0,p2) > 
da(0,Pi) + c/2. 

Recall that O and P2 are two points on the surface of a ball B in the metric 
d with centre Q2 and radius r2 = Q(k 1 / 3 ). We claim that the midpoint M of 
OP2 is 'well inside' this ball, i.e., that 

d{M, Q 2 ) < (1 - s)r 2 (4) 

for some constant e > depending only on d, A and C . Recall that d(0, P2) > 
ak 1 / 3 and r 2 < ^4fc 1/3 , so d(0,P 2 ) = 6(r 2 ). If d = d 2 , then (gj) follows imme- 
diately: in fact, it follows immediately for any metric d defined by a 'strictly 
convex' norm, i.e., one whose unit sphere contains no line segments. 

For d = c?jm, (H) follows using the additional fact that the angle between 
OP2 and the plane is not within c\ of 7r/4: the unit sphere of the corresponding 
norm is a double cone. The only line segments this contains are at an angle 
7r/4 to the plane. Let H be the set of pairs (A, B) of points on this unit sphere 
with d2{A : B) at least some small constant and the angle between AB and the 
plane lying outside the interval (7r/4 — c\, 7r/4 + ci). Then \\{A + P)/2||jm < 1 
for all (A, B) e H, and hence, by compactness, \\(A + P)/2||jm < 1 — £ for all 
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Figure 3: Two points Pi, P2 of X such that the angle P\OP 2 is very small. P[ is 
the point of OP2 at distance d(0,P\) from O. As P1OP2 is small, the distance 
PiP[ is o(l). Since P 1 P 2 > c, it follows that P[P 2 > c/2. The points and P 2 
are equidistant from Q2; the dotted line is part of the sphere d{Q 2 , •) = d{Q 2l O) 
(curvature exaggerated), and M is the midpoint of OP 2 . Using convexity, we 
can show that P\ lies inside the corresponding ball. 

(A, B) 6 H, for some e > 0. Relation g]) follows, taking A = (O - Q 2 )/r 2 and 
B=(P 2 -Q 2 )/r 2 . 

The point P[ lies on the line segment OP 2 , and (recalling that d 2 (0 : P\) > 
afc 1 / 3 ) is at Euclidean distance at least c/2 from both endpoints of OP 2 . Thus, 
P{ lies on XM where X = O or X = P 2 , and d 2 (X,P{) > c/2, so A = 
g?2(X P{)/d 2 (X, M) > C2fc -1 ^ 3 for some constant C2. From convexity of d we 
have 

d(P{,Q 2 ) - d((l-A)X + AM,Q 2 ) 

< (l-X)d(X,Q 2 ) + Xd(M,Q 2 ) 

< (1 - A)r 2 + A(l - e)r 2 =r 2 - e\r 2 = r 2 - 9(1). 

Since d(Pi,P[) = o(l), it follows if k is large enough that d{P\,Q 2 ) < r 2 , 
contradicting our assumption that d(Pi, Qj) > fj for all i and j. □ 

The proof above can be easily adapted to any norm on M 3 by redefining the 
event -B4: we must exclude a set of directions that form a neighbourhood of the 
set of directions in which the unit sphere contains a line segment. One can check 
that this latter set has measure zero for any norm, so it has neighbourhoods 
with arbitrarily small measure, and we can ensure that P(i?4) = o(e~ ck ) as 
above. 

We now turn to the proof of Theorem [8l The key observation is that a 
point of the s by s torus T(s) is robustly black with respect to (P 2 ,V 2 ) if 
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it remains black when the points of V£ are shifted away from the plane by a 
suitable distance. As a shifted Poisson process is again a Poisson process (but 
on a different set), this allows us to construct the coupling so that 'defects' arise 
essentially independently - this is the reason why Theorem [5] is easier to prove 
than the corresponding result for random Voronoi tessellations in the plane. 
Theorem 6.1 of [3]. More precisely, we can deduce Theorem [5] from the ungainly 
Lemma 1101 below. In this result, the notation is rather unnatural - we write 
V\ and H.2 for coloured Poisson processes where the probability that a point 
is black is pi and P2, respectively. We write IZ2 rather than V 2 since we shall 
modify the process IZ2 to obtain a process V2 with the properties required for 
Theorem |U 

Lemma 10. Let < pi < P2 < 1 and e' > be given. Set §' = s~ £ . We may 

construct in the same probability space Poisson processes , , and TZ^ on 
T(s) x [0, s] of intensities pi, 1 — pi, and 1 — p 2 , respectively, and a Poisson 
process 1Z\ on T(s) x [5' , s] of intensity P2, so that "Py and are independent, 
IZ2 and 7Z2 are independent, and the following global event E' gl holds whp as 
s — > 00: for every piecewise-linear path Pi C T(s) which is black with respect to 
CPi^i) there is a piecewise-linear path P2 C T(s) which is black with respect 
to (Jii, , IZ2 ) , such that every point of P2 is within distance log s of some point 
of Pi and vice versa. 

Before proving Lemma[l0l let us show that it implies Theorem[8l From here 
on we follow the convention in [3] of writing, x, Xi etc for points of T(s), and 
z, z' etc for points of our Poisson processes in T(s) x [0, s]. (The more natural 
notation P G V used above becomes confusing when paths P are involved.) 
With this convention a point x is black if the nearest z is black. 

Proof of Theorem\^ We may assume without loss of generality that 6 = s~ e < 
1. Set e' = e/3, and note that 8' = s -e / 3 is larger than S. Let Vi , "P-f , TZ J and 
IZ2 be coupled Poisson process with the properties described in Lemma [TUl set 
T2 = 1^2 ' an d l et 

Vt = {(x, y, z-S'): (x, y, z) e TZ+} U V , 

where V' is a Poisson process of intensity P2 on T(s) x [s — S',s] that is inde- 
pendent of Tlf. Note that and V2 are independent Poisson processes on 
T(s) x [0, s] with the desired intensities p 2 and 1 — p 2 - It remains only to show 
that E g \ holds whp. 

From Lemma fTOl we may assume that E' gl holds. Let E be the event that 
every point of T(s) is within distance ©((logs) 1 / 3 ) = 0*(1) of some point of 
1^2 , say. (As usual, we write f(s) = 0*(g(s)) if f(s) = 0((log s) c g(s)) for some 
constant C.) Then, from basic properties of Poisson processes, E holds whp. 
To complete the proof, we shall show that E' gl (~l E implies E s \. 

Suppose that E' gl and E hold, and let Pi be a path that is black with respect 
to {Vt,Vi ). Then, since E' gl holds, there is a path P2 within Hausdorff distance 
logs of Pi that is black with respect to (IZ^TZ^)- Let x be any point of P2, 
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and let z + and z~ be points of TZ% and TZ^ at minimal distance from x. Since x 
is black and E holds, we have d(x, z + ) < d(x, z~) = 0*(1). Let z' be the point 
of V2 obtained by shifting z + by a distance 6' in the negative z-direction. We 
claim that 

d(x,z') < d(x,z + ) - (S') 2 /0*{1). (5) 

If d = djM, this is immediate; the reduction in distance is exactly 6'. If d = cfei 
then the extreme case is when z 1 lies in T(s) and x and z' are at maximal 
distance, in which case the claim follows from Pythagoras' Theorem. 

Now (<5') 2 = s~ 2e / 3 , while 6 = s~ 6 . If s is large enough, then as = TZ^ , 
it follows from © that 

d(x,V£) < d(x, z') < d(x, z+) - 2C d S < d(x, z~) - 2C d S = d{x,V^) - 2C d 6. 

Thus x is (2C£i(5)-robustly black with respect to (V^iV^)- Since x was an 
arbitrary point of P2, it follows that P2 is (2Cd^)-robustly black. Finally, as Pi 
was arbitrary, E g \ holds, as required. □ 

We now turn to the proof of Lemma [TOl 

Proof of Lemma \l(A Throughout the proof we write 5 and e for the quantities 
8' and e' appearing in the statement of the lemma. 

To construct our coupled Poisson processes, we start with three independent 
Poisson processes, a process V of intensity 1 on T(s) x [5, s], and two processes 
on T(s) x [0, 5], a process Vs of intensity 1 and a process V of (much higher) in- 
tensity J" 1 / 2 . We shall form (homogeneous) Poisson processes Vi and TZf with 
the properties described in the statement of the lemma by assigning every point 
of V to exactly one of V^ 1 and to exactly one of IZ^ , assigning certain points of 
Vs to Vi and/or TZ^ , and assigning certain points of V to Vf ■ [Thinking of 
V~i as a single Poisson process whose points are coloured black/ white according 
to a colouring col 1; and writing TZ^ similarly in terms of (7^-2, cohj), then every 
point z of V is present in both V\ and 7?-2, although the colours coli(z) and 
coi2(z) may be different. Certain points of Vs are present as white points in V\ 
and/or TZi, and certain points of V are black in V^ .] 

Table [T] summarizes the domains and intensities of these Poisson processes, 
as well as two others that we shall consider in the proof. The first four lines 
show the processes we shall construct; the remaining lines concern processes 
used in the construction. 

Table [2] shows the probabilities with which the points of our independent 
processes V, Vs and V are included into the derived processes. Note that V^ 
is simply V^ s _ s U V. 

In constructing V^ and TZ^ we shall define two intermediate processes, 
Vi s _ s and T>, whose union will form V^ . The notation reflects the fact that 
V{~ s will consist of almost all of V^ , more precisely, all points except those 
with z-coordinate at most 5. We write T> for the remaining points of V^ since 
these points will form 'defects' in our first attempt at a coupling. 
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Process 


Domain 


Intensity 




T(s) x [0, s] 


Pi 




T(s) x [0,3] 


1-P1 




T(s) x [6, s] 


P2 




T(s) x [0, s] 


1 ~ P2 


V 


T(s) x [<5, s] 


1 


Vs 


T(s) x [0, S] 


1 


V 


T(s) x [0,(5] 




' l,s-S 


T(s) x [5, s] 


Pi 


V 


T(s) x [0,(5] 


Pi 



Table 1: The various Poisson processes involved in the statement and proof of 
Lemma [TU1 We must construct and to be independent, and 7Z 2 and TZ^ 
to be independent. The processes V, Vs and T>' are independent by definition; 
Vi s _ s and T'-f will be independent by construction. 





V 


v + s 

1 l.s-6 


vt 


k n 


1Z 2 


V 




Pi 


Pi 


1 - Pl P2 


1 -P2 


Vs 








1-Pl 


1 -P2 


V 


PiS 1 / 2 




PK5 1 / 2 







Table 2: The matrix of inclusion probabilities when X>, V^ s _ s , V^ 1 and 72.^ are 
constructed from V, Vs and V. Note that 7?+ = V^ s _ s U P. 

By the 'natural' coupling of V 1 with TZf we shall mean the coupling obtained 
as follows. Given V, Vs and V as above, for each point z of V, toss a three- 
sided coin: with probability pi assign z to both V{~ g and 7Z 2 , with probability 
P2 — Pi assign z to 1Z 2 and T'] - , and with probability 1 — p 2 assign z to 7Z 2 and 
Vi . Similarly, for each point z of Vs, with probability p\ assign z to neither 
Vi nor 72.^" , with probability p 2 — pi assign z to P-j - , and with probability 
1 — pi assign z to both P-f and 1Z 2 . Given V and Vs, we make all these 
choices independently. Thus, the sets 1Z 2 and 1Z 2 we obtain are independent 
Poisson processes on T(s) x [S, s] and T(s) x [0, s], respectively, with respective 
intensities p 2 and l—p 2 . Also, V^~ s _ s andVi are independent Poisson processes 
on T(s) x [(5, s] and T(s) x [0, s] respectively, with respective intensities pi and 
1 — pi. To obtain V^ , form a set T> by selecting each point of V independently 
with probability piS 1 / 2 , and set Vf = s U V. 

The construction above has the property that Vi a _ s C 1Z 2 and D 1Z 2 . 
Thus, any path Pi in T(s) that is black with respect to (Vi s _ s ,Vi) is black 
with respect to (1Z 2 , 1Z 2 ) . Unfortunately, when we add the points of T> to Vi s _ s 
to obtain T^, this may introduce new black paths with respect to (V^,Vi), 
which need not be black with respect to (1Z 2 , TZ 2 ) . For this reason, we think of 
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the points of V as defects. To deal with these defects, we shall need to adjust 
the coupling. 

By a potential defect we mean a point of T>' . Each potential defect has only 
a small probability, piS 1 ^ 2 , of becoming a real defect. Even though there are 
many more potential defects than defects, the density of potential defects is still 
low: the projection of V onto T(s) is a two-dimensional Poisson process with 
intensity S^^S = 5 1 ' 2 = s~ e l 2 . 

In constructing our final coupling, we shall condition on T 5 , Vg and £>'. Given 
these processes, we shall construct the remaining processes in a way that respects 
the inclusion probabilities shown in Table [U the only changes we shall make 
to the 'natural' coupling just defined are to the coupling of inclusion choices 
associated to with inclusion choices associated to Tlf . The first step is to 
describe certain very unlikely 'bad' events defined in terms of V, T$ and V . 
When one of these events holds, we shall complete the coupling arbitrarily (for 
example, as above), and E' gl will not necessarily hold. As the bad events will 
have probability o(l), this will not be a problem. 

Let A be a large constant, and let B\ be the event that there is some point of 
T(s) for which no point of V lies within ci-distance A(logs) 1 / 3 . An elementary 
calculation shows that P(i?i) = o(l) if A is chosen large enough. (For d = di, 
we may take A = 1, or indeed any constant such that ttA 3 /3 > 1.) 

Let us say that two potential defects, i.e., points of V , are close if they are 
at d-distance at most 4A(logs) 1 / 3 ; this relation may be taken to define a graph 
on V . By a cluster of potential defects we mean a component of the resulting 
graph. Let B2 be the event that there is a cluster of potential defects containing 
more than 10/e potential defects. If B2 holds then, considering a connected 
set of exactly [10/e] potential defects, there is a d-disk D C T(s) of radius 
r = [10/e]4A(logs) 1 / 3 = 0((log s) 1 / 3 ) = 0*(1) containing the projections of at 
least [10/e] potential defects. But we may cover T(s) with 0(s 2 ) <i-disks D L of 
radius 2r so that any d-disk of radius r is contained in some Di. The projection 
of V onto T(s) is a Poisson process of intensity S 1 ^ 2 = s~ e / 2 , so the expected 
number of images in a given Di is A = S 1 ^ 2 area(D i ), which is at most s~ £ / 3 if 
s is large enough. Hence, 

¥(B 2 ) < P(3A containing > [10/e] images) 

* °( s2 )^! e "- O ( s2Ari0/el )-°( 1 )- 

Let us say that a point z' of T>' and a point z of V are potentially adjacent 
if the Voronoi cells of z and z' in the tessellation of T(s) associated to the point 
set V U {z'} and metric d meet. Note for later that whenever the Voronoi cells 
of z G V and z' € V defined with respect to some point set X containing 
V U {z 1 } meet, then z and z' are potentially adjacent: deleting points of X to 
obtain V U {z 1 } can only enlarge the Voronoi cells associated to z and z'. We 
shall apply this observation later with X = T> U V C U V\ ■ Note also that 
whether or not z and z are potentially adjacent depends only on "P, V5 and 2?', 
not on Pj^, which we have not yet constructed. 



19 



Let a be a small constant to be chosen later, and let B3 be the event that 
there is some z' £ V that is potentially adjacent to at least a logs points 

Z!,Z 2 ,...,Zk of V. 

Claim. For any choice of the constant a > 0, we have ¥(B^) = o(l) as s — > 00. 

Since P(Si) = o(l), the claim follows if we show that P(BJ n B 3 ) = o(l). 
Wc shall deduce this from Theorem [9] 

For z' £ 2?' , let -E^/ be the event that there are k = \a log s] points Zi, . . . , Zk 
of P and corresponding points Xi, . . . ,Xk £ T(s) x {0} such that d(z',Xi) = 
d(zi,Xi) = Ti < A(logs) 1 / 3 , and d(zj,Xi) > ri for all i, j. If B\ n -B3 holds, 
then so does i? z / for some z' GT>': there is some z' £ V potentially adjacent to 
k points Z\,...,Zk £ V. The Voronoi cells of z' and Zi defined with respect to 

V U {z'} meet at a point xi that is equidistant from z' and Zi, with 2) > 
d(xi, Zi) = d(xi, z') for every point z £ V, and in particular for all Zj. Since B\ 
holds, some point of V is within distance A(logs) 1 / 3 of Xi, so r*j = d(xi,Zi) < 
A(logs) 1 / 3 , and .B^ holds. 

Note for later that, if £?i does not hold, any two potentially adjacent points 
are within distance 2A(\og s) 1 / 3 . 

To show that P(Bl n S 3 ) = o(l), we shall condition on V. Note that 
E(|X>'|) = S~ 1/2 Ss 2 < s 2 /2. Since \V'\ has a Poisson distribution, it follows that 
\D'\ < s 2 whp; in proving the claim we may thus assume that V is fixed and 
that \V\ < s 2 . Let z' be a point of V . As V and 2?' are independent, extending 

V to a Poisson process on T(s) x K and translating through the vector —z', so 
that z' is moved to the origin, we may realize V ' ■ X 0jS 0/ subset of a Poisson 
process on T(s) x M of intensity 1. In this Poisson process, if E z > holds, then 
the shifted points Pi = Zi — z' and Qi = Xi — z' have the properties described in 
Theorem [SJ with k = [a logs]. Of course, this result concerns a process on R 3 
rather than on T(s) x R. But the event considered only involves points within 
distance (^((logs) 1 / 3 ) = o(s) of the origin, so this makes no difference. Hence, 
by Theorem H P(E Z ,) = o(e" 3fc / a ) = o(s~ 3 ). Considering the 0(s 2 ) points of 

V separately, it follows that ¥(B\ n B 3 ) = o(l). As P(Bj) = o(l), this proves 
the claim. 

Let us say that a point z' £ V and a point z £ Vg are very close if they 
are within distance 2A(logs) 1//3 . Let B4 be the event that some z' £ V is very 
close to at least a logs points z £ Vg. As the projection of Vg onto T(s) has 
intensity 5 = s~ e , it is easy to check that whp no point of V is very close to 
more than 10/e points ofVg. (The argument is similar to but simpler than that 
for P(£?2) = o(l) above, so we omit the details.) It follows that P(£?4) = o(l). 

From now on we condition on V, Vg, and X>'; we regard these sets as fixed 
for the rest of the proof, and assume, as we may, that none of B±, . . . , B4 holds. 
To complete the construction of the coupling, it remains to assign each point z 
of V to V^ s _ s with probability pi and otherwise to Pf, to assign z to TZ^ or 
1Z.2 with probabilities P2 and 1 — p 2 , to assign each point z of Vg to with 
probability 1 — p\ and to TZ^ with probability 1 — to select points of V with 
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probability pid 1 ^ 2 to form V, and then to set V^ = V^ 8 _ s U V. As long as 
all choices involved in constructing (Vi 7 V^~) are independent, and V^ will 
have the right marginal distribution; the same holds for (TZ^ , TZ^). ^ we can 
also ensure that E' gl always holds, the coupling will have the properties claimed 
in the statement of the lemma. 

Let Ci, C2, . . . C V be the clusters of potential defects. By assumption, B2 
does not hold, so no Ci contains more than 10/e points of V . Let r(Ci) denote 
the set of points z of V U Vs with the property that z is potentially adjacent (if 
z G V) or very close (if z G Vs) to one or more points z' G Ci. As B\ does not 
hold, any pair of potentially adjacent points is at distance at most 2A(\ogs) 1 ^ 3 . 
Any pair of very close points is also separated by at most this distance. Hence, if 
z G T(Ci)nT(Cj), then the clusters Ci and Cj contain close points, contradicting 
the definition of a cluster. Thus, the sets r(Ci) are disjoint. Also, as B3 U B4 
does not hold, we have |r(C^ ) | < 2a\ogs\Ci\ < 20ae~ 1 logs for every i. 

We shall construct our coupling independently for each set Ci L) T(Ci) in 
a manner described below. For each point z £ P not in [Ji^(Ci), we toss a 
three-sided coin, assigning z to V^ s _ s and 1Z\ with probability p\ , to 'P-f and 
IZ2 with probability P2 — pi, and to V\ and TZ^ with probability 1 — pi . These 
choices are independent for different z, and also independent of the choices made 
within the Ci U T(Ci); this corresponds to the 'natural' coupling described at 
the start of the proof. Similarly, for each z G Vs not in [J 4 T(Ci), we include z 
into both Vi and TZ^ with probability 1 — p-2, and into Vi only with probability 
P2 - Pi- 

Let C be a cluster of potential defects. Let V, V^ s _ s , Vf~ and IZf be 
defined as in the natural coupling, but restricting our attention to C U r(C), 
i.e., starting from C, T(C) n Vs and T(C) D V in place of V, Vs and V. Note 
that V^ and TZ^ have the correct distributions for the restrictions of Vf 1 and 
K^toCU F(C); however, to ensure that E' gl holds we shall have to adjust the 
coupling while keeping the marginal distributions fixed. 

Let B(C) be the 'bad' event that T> is non-empty, i.e., that one of the 
potential defects in C is an actual defect. Since |C| < 10/e, and each z' G C is 
included in T> with probability pid 1 / 2 = p\s~ e / 2 , we have P(B(C)) < s~ £ / 3 for 
s large. 

Let G(C) be the 'good' event that every point of T(C) is in V^[ but no point 
of T(C) is in TZ~2. (Thus, each z G Y{C) n V is in ^J, while each z G r(C) n Vs 
does not appear in TZf ■) From the definition of the natural coupling, G(C) 
has probability (p 2 -pi) |r ' c) '. Recall that |r(C)| < 20ae~ 1 logs. By choosing 
the constant a sufficiently small, we may ensure that (p2 — pi)' r ^ c '' > 2s~ £//3 . 
Hence, P(G(C)) > 2P(B(C)), and there is some event G'(C) G G(C) \ B(C) 
with probability P(B(C)). [To be pedantic, we must modify our probability 
space at this point. Having conditioned on V, Vs and V, and restricting our 
attention to coordinates involving points in CL)T(C), we are working in a finite 
probability space. To ensure that we can choose an event G'(C) with exactly 
the right probability, we should work in a space without atoms, so we simply 
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adjoin one extra random variable Uc for each cluster C, with Uc uniform on 
[0,1], say, and the new variables independent of everything else.] 

Wc define our final coupling by 'crossing over' the natural coupling on the 
events B(C) and G'(C): writing Q for the probability space on which we have 
defined and TZf , let / be a measure preserving bijection from B(C) UG'(C) 
to itself, mapping B(C) into G'(C) and vice versa. We define the restrictions 
of to C ur(C) to be equal to . Temporarily abusing notation by writing 
Tlf for the restriction of K± to CL)T(C), for cu <£ B(C)UG'(C) we set nf(to) = 

TlJ (w), while for cu E B(C)UG'(C) wc set IZf(cu) = nf{f(u;)). We extend the 
coupling to all clusters independently. As the crossing over does not affect the 
marginal distribution of (the restriction to C U T(C) of) 7tL^, we end up with 
the correct marginal distributions for V^ 1 and Tlf ■ It remains to check that E' gl 
holds. 

We claim that in our final coupling 

T>t,a-S C n 2 and V X D n 2 ( 6 ) 

always holds. The 'natural' coupling has this property, so we must check that 
it is preserved by the local 'crossing over' within C U T(C). Recall that C C V 
consists only of potential defects, which never appear in any of the processes in 
([6]): some potential defects will be selected to form V, which will later be added 
to Vi a _ s to form Vf. Recall also that L(C) C V U Vs. Finally, recall that we 
only 'cross over' at elements w £ B(C) U G'(C) of the state space. 

Suppose first that u E G'(C) C G(C) \ B(C). Then, since w E G(C), every 
point of r(C) is in , and (hence) none is in s _ s . Thus the restriction of 
gtoCU r(G) holds in this case. 

Suppose next that w E B(C). Then f(u) E G'(C) C G(C), so every point 
of r(C)fl? is in 72-2 and no point of T(C) is in TZ^ ■ Recalling that points of Vg 
never appear in Pi s _ s , it again follows that the restriction of © to C U T(C) 
holds. Since C was an arbitrary cluster, this establishes ([6]). 

From ([6]), if x E T(s) is black with respect to (V^~ , 'P-f ) and the closest point 
of Vt U T>i is not a defect (a point of V = \ V^ s _ s ), then x is black with 
respect to (TZ% , TZ^ ) ■ It remains to show that we can modify a black path to 
avoid the defects; we shall do this by taking a short detour around each defect. 

For each defect z', let V z > be the Voronoi cell of z' in the tessellation V of 
T(s) defined with respect to T> U V C Vi U . Note that V z i contains the cell 
of z' in the tessellation associated to (V^,Vi). Let {z[, . . . , z' r } C V C V be 
a maximal set of defects such that the union U of the cells V z > is connected, 
and let d°°U denote the external boundary of U. As B\ does not hold, no 
cell V z i has radius larger than A(logs) 1 / 3 , so the z[ lie in a single cluster C of 
potential defects. Since B^ does not hold, it follows that U has diameter at 
most ©((logs) 1 / 3 ). 

Claim. Every point of d°°U is black with respect to (7^.^,7?.^). 

Before proving this claim, let us note that it implies that E' ^ holds. Let P\ be 
a piecewise-linear path in T(s) which is black with respect to (P^ , P-f ) . Every 
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point of Pi not in one of the sets U is black with respect to (TZ^ , TZi ) > since 
the nearest point of V± U Pf is not a defect. On the other hand, every point 
of some set U is surrounded by a black cycle, d°°U , that lies within distance 
logs. Thus, whenever Pi visits a set U, we may replace a part of Pi by a part 
of d°°U , obtaining a path Pi that is black with respect to TZ^), with every 
point of Pi within distance log s of some point of P2 , and vice versa. Shifting 
Pi slightly, we may assume that Pi is piecewise linear. As Pi was arbitrary, 
this shows that E' gl holds. 

All that remains is to prove the claim. Let x G d°°U. Then x is in some 
cell V z of V with z £ {z[, . . . , z' r }. Since U U V z is connected, the maximality 
of {z[, . . . , z' r } implies that the point z is not a defect. Therefore, recalling that 
the tessellation V is defined with respect to V UP, we have z G P. Furthermore, 
no point of P is closer to x than z is. In the tessellation V, the cell Vz meets one 
of the cells V z >, z[ G C. Hence, z is potentially adjacent to z[, so z G r(C). As 
the point z{ G C is a defect, B(C) holds, i.e., w £ S(C). Hence, /(u>) G G(C). 
Since z G r(C) (~1 P, it follows that z G 72^. It remains to check that no point 
z" of IZi is closer to x than z is. As 7?.^ C P U Pj, and z is the closest 
point of P to x, it suffices to consider points z" G Vs. But 73i holds, so 
d{x,z$ = d(x,z) = d(x,V) < A(logs) 1 / 3 , so any z" G P 5 withd(cc,z") < d(cc,z) 
is very close to z[, and so lies in T(C) and hence, from the definition of G(C), 
cannot be in 1Z% ■ Hence, every x G d°°U is indeed black with respect to 
(7Zi,lZi), completing the proof of the claim, and hence of Lemma 1X171 □ 

5 Extensions 

For simplicity, Theorem [T] was stated and proved for tilings of K 2 generated 
by two specific norms on R 3 , namely djM and di. In the context of growing 
crystals, it is natural to consider certain other norms, for example the ^i-norm 
on R 3 , corresponding to the sum of £1 on R 2 and time. The statement and 
proof of Theorem [1] adapt immediately to this setting: the only change required 
is the definition of the event -B4 in the proof of Theorem [9l In fact, the same 
comment applies to any norm of the form 1 1 • 1 1 ffi t\, where 1 1 • 1 1 is a norm on R 2 
which has the symmetries of the square; these symmetries are assumed when 
we apply the Russo-Seymour- Welsh type result from [3] 

Theorem[T]and its proof also extend to two-dimensional slices of d-dimensional 
(usual) Voronoi tessellations, this time with simple changes: for example, the 
bound ©((logs) 1 / 3 ) on the distance to the nearest point of P must be replaced 
by 0((log s) l ' d ), and T(s) x [0, s] is replaced by T(s) x [0, s] d ~ 2 . When we 'shift' 
points, it suffices to change one coordinate, so the set T(s) x [6', s] appearing in 
the statement of Lemma fTTJI is replaced by T(s) x [5' , s] x [0, s] d ~ 3 . 
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